Molecular orientation entanglement and temporal Bell-type inequalities 
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We detail and extend the results of [Milman et al. , Phys. Rev. Lett. 99, 130405 (2007)] 
on Bell-type inequalities based on correlations between measurements of continuous observables 
performed on trapped molecular systems. We show that, in general, when an observable has a 
continuous spectrum which is bounded, one is able to construct non-locality tests sharing common 
properties with those for two-level systems. The specific observable studied here is molecular spatial 
orientation, and it can be experimentally measured for single molecules, as required in our protocol. 
We also provide some useful general properties of the derived inequalities and study their robustness 
to noise. Finally, we detail possible experimental scenarii and analyse the role played by different 
experimental parameters. 

PACS numbers: 03.65.Ud;03.67.-a;33.20.Sn 



I. INTRODUCTION 



Quantum mechanics allows for the existence of states 
without any classical correspondence. Examples of such 
states are entangled states, that may appear when de- 
scribing the total state of a many particle system, or 
when describing different degrees of freedom of a sin- 
gle particle. Entangled states cause debate because they 
can present properties which contradict our classical intu- 
ition. In addition, some of these properties can be use to 
increase the security of quantum communication and the 
efficiency of algorithmic protocols when compared with 
classical techniques. This is why so much attention has 
been payed to problems belonging to the foundations of 
quantum mechanics, and in particular to quantum entan- 
glement, which is believed to be one of its main traits 0. 
One of the most important discussions on entanglement 
and its conflicts with classical physics concerns the real- 
ism and (non-)locality of quantum physics. As pointed 
out by Einstein, Podolski and Rosen in the so called EPR 
[H paradox, entangled states are closely connected to ap- 
parent contradictions between quantum mechanics and 
fundamental physical assumptions. Even if experimen- 
tal evidence has been obtained to support quantum me- 
chanics against the EPR criticism it is still a matter 
of debate if such experiments were realised in the ideal 
conditions, closing all the loopholes, so as a definitive 
conclusion can be reached. At the same time, from the 
fundamental point of view, identifying precisely which es- 
sential quantum properties are involved in local realism 
violation is still an open problem. This is why extend- 
ing local realism tests to different physical systems and 
different physical scenarii still presents so much interest. 

A common property of non-local states is entangle- 
ment. In spite of its importance and consequences in 
different physical contexts, it remains an open question 
to determine whether a general quantum system is en- 
tangled or not and to quantify the degree if entangle- 
ment of a given state. The problem has been solved for 
some particular cases, as for a general bipartite system 



of dimension H2 ^ ^2 and '^^2 ® '^3 , where Hd is the 
one particle Hilbert space of dimension d. In such cases, 
necessary and sufficient conditions for telling whether a 
given state is entangled or not exists. In other particular 
cases, or subspaces, one can also find necessary and suffi- 
cient conditions. For instance, entanglement in bipartite 
pure states can always be recognised and quantified ir- 
respectively of each parties' dimension. However, when 
dealing with arbitrary states, including the more realis- 
tic mixed ones, only necessary conditions for separability 
(non-entanglement) can be provided. A notion that will 
be useful in the following of this paper is the one of en- 
tanglement witnesses, defined as an operator W for which 
the expectation value {W) < S for all separable states. 
This ensures that the state is entangled if (W) > S. On 
the contrary, the case {W) < S ^ is not conclusive. Ex- 
amples of entanglement witnesses that are also useful for 
fundamental tests of quantum mechanics are Bell-type 
inequalities (l^ , which are the main scope of this paper. 

Bell-type inequalities are composed of combinations 
of observables that, when measured, allow for setting a 
board between crucial aspects of quantum theory and the 
classical one. They were formulated by J. S. Bell 0] as a 
late reply to the EPR criticism to quantum physics 
Bell inequalities aim at answering the question: is quan- 
tum mechanics a local and realistic theory? In order 
to do so, they combine correlations between local mea- 
surements realised in a multiparty system. Since their 
first formulation, several other inequalities have been pro- 
posed, some studying the same type of problem as Bell, 
other generalising locality and realism tests to many ob- 
servers and many possible observables. There are a num- 
ber of Bell type inequalities, and classifying all of them 
is a work by itself [8| . 

The inequalities derived here follow the original for- 
mulation of Clauser, Hornc, Shimony and Holt (CHSH) 
[J, and deal with the following scenario: two observers A 
and B perform local measurements on a bipartite system. 
Each observer can chose among two experimental set-ups 
(a and a' for A and b and b' for B). Each measurement 
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performed by A and B can give only two outcomes. At 
its origin, CHSH inequalities have been formulated for 
a pair of spin 1/2 particles or equivalent two-level sys- 
tems. In the framework of local hidden variable (LHV) 
theories the measurement outcomes correlation statistics 
which must fulfil: 

\{o-aO-b) + {o-aCTb') + (c^a'Crb) - {cra'Crb') \ < 2, (1) 

where ctq is the Pauli matrix in the a direction. Briefly, 
in a LHV theory, one assumes that measurements per- 
formed by each observer are independent and their out- 
comes have a probability distribution which is a product 
of independent probabilities for each subsystem. Such 
probabilities can also depend on some random local vari- 
able. Details are discussed in many works, as [2^, for 
instance. It can be shown that some entangled states 
can violate ([1]), and this experimental violation was ob- 
served using photon pairs entangled in polarization Q. 
In this case, directions a, a', b and b' refer to different 
orientations of polarizers placed before the detectors, de- 
termining the direction of the Pauli matrix that is mea- 
sured. It was shown by Cirelson (29j that the maximum 
value of ((T|) is 2 a/2, and it is easy to verify that this max- 
imal violation can be obtained with maximally entangled 
states. 

Inequalities as ^ have proved to hold for two-level 
systems or equivalent ones. By "equivalent ones" , we 
also include Bell type inequalities involving a continuum 
of possible measurement outcomes that are dichotomised 
and transformed into a two outcomes measurement set- 
up. Dichotomization works as follows: one splits in two 
classes the range of possible measurement outcomes. All 
results obtained lying in one of the classes is identified to 
a given value (+ or — ) and the results obtained in the 
complementary space are associated to the opposite sign. 
Some examples of systems where this can be sucesfully 
done are optical fields 0, . An interesting problem is 
to derive Bell-type inequalities for continuous variables 
without appealing to dichotomisation. Cavalcanti et al. 
found a way out by using second moment correlations 
instead of first moment ones, as done in CHSH- type in- 
equalities as (1) [lO]. Here, we deal with this problem 
in a different way: by using bounded observables, one 
can still use CHSH-type inequalities to detect non-local 
properties. In this Cirelson bound depending on 

the norm of the measured observable can also be derived, 
even if, at least for the specific case treated in this pa- 
per, we have not shown yet that it can be attained. An 
interesting property of the inequality discussed in this 
work is that it can be used not only for infinite dimen- 
sional systems, but also in N levels systems, where N 
is a finite number. In this case, the numerical value of 
the bound splitting between a local theory and non-local 
one depends on the maximal eigenvalue of the measured 
observable. 

Up to now, several studies have been made on Bell-type 
inequalities in different contexts. The general conclu- 
sion is that the subject still presents several open ques- 



tions and no general theory is available. In particular, 
a number of intriguing features coming out from such 
studies somewhat contradict our acquired quantum "in- 
tuition" : in [2H it is numerically shown that bipartite 
multidimensional states may violate locality tests more 
than two qubits. For two qubits, it has been proven 
that a maximal violation exists, and it is given by the 
Cirelson bound [2^, as will be discussed hereafter. Acin, 
Gill and Gisin showed, some years later, also using nu- 
merical tools, that the maximal violation for pure bipar- 
tite multidimensional systems is not obtained for maxi- 
mally entangled sates [2J|. General rules relating local 
realism violation and entanglement have not yet been 
found, even if it has been proven that all non-local states 
are entangled in some way. This fact can be easily un- 
derstood with the help of entanglement witnesses. It has 
been shown in [12| that Bell-type inequalities are entan- 
glement witnesses, and all non-local states are entangled. 

The inequalities studied in this paper are based on 
molecular spatial orientation correlations measurements 
[2^ instead of spin-like observables. However, the same 
type of idea can be generalised to other continuous 
bounded observable. We have shown that they can be 
implemented using time delayed measurements of corre- 
lations between the spatial orientations of two molecules. 
As in usual Bell inequalities scenarios, the proposed non- 
locality tests rely on measurements performed indepen- 
dently on each molecule by observers placed far apart 
enough, so that no communication between them is possi- 
ble during the realization of the protocol. We have shown 
that the proposed inequalities can be violated by a set 
of entangled states. An interesting point is that the in- 
equalities derived in the present paper can also be used 
as entanglement witnesses, as shown, for general CHSH- 
tyc inequalities in Q. In this situation, one can loosen 
measurement conditions, since we are interested only in 
detecting a particular quantum correlation, and not a 
fundamental aspect of quantum physics. 

From the experimental point of view, motivations for 
this work are the recent advances in single molecule ma- 
nipulation and detection with entanglement creation [Tsj 
and quantum information purposes [13, [H, [O, [ill . In 
particular, trapped cold polar molecules are promissing 
candidates for quantum information processing based on 
the manipulation of their rotational levels [23]. Rota- 
tional states, which are also involved in the Bell-type 
inequalities studied in the present paper, are relatively 
long lived, allowing for short quantum gate implemen- 
tation times: one can hope to perform about 10* gate 
operations before decoherence takes place. This excel- 
lent performance when compared to cold collision based 
quantum gates [l^ is due to the strength of molecular 
interactions, based on dipolar forces. 

The present paper is organised as follows: in section HIl 
we discuss some general properties of Bell-type inequal- 
ities that are useful in the context of the non-locality 
tests we propose. In section Hill we explicit the inequali- 
ties based on molecular orientation and show how it can 
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be used for non-locality tests. We exploit its performance 
and study some entangled states that violate it. We then 
discuss some constraints involved in an eventual physi- 
cal implementation of our ideas in section [TVl ending up 
with a concluding note in section FVl 

II. GENERAL PROPERTIES 

We describe now some general properties of the in- 
equalities studied in the present paper. Our inequalities 
involve four observables that can be combined as follows: 

B= Ol(0l)(»O2(02)+Ol(</'l)®O2(0^)+ (2) 

Oi(0'i) ® 02(02) " Oi(0;) <^ 02(0^). 

The operators appearing in the equation above are de- 
fined as follows: operators Oi{(j)i), i = 1,2 are local ob- 
servables chosen by observers 1 and 2. The variables 
(pi also depend on local properties only, and Oi{4>i) = 
U{(j)i)0{Q)U\(j)i), with U{4>i) a one-parameter group of 
unitary operators with periodicity 27r. If the operator 
Oi{(j)i) is bounded, we can show that, under the assump- 
tion of local realism, ([2]) satisfies: 

m\<s, (3) 

where S* is a number, representing the maximal allowed 
value for a local theory to hold whenever measurements of 
correlations between observables are compared as in ([2]). 
The basic assumption to derive Q is that correlations 
can be described by probability distributions which are 
independent for each party (1 and 2) and depend only on 
local parameters and that the spectrum of observables 
Oi[(j)i) is bounded. The same inequality can be derived 
by assuming that the average of B is taken with respect 
to a separable (non-entangled) states. 

For simplifying reasons, we focus on the specific case in 
which the norm of Oi{(f)i) is the same for each subsystem 
and is equal to Xmax- In this case, the numerical value 
of5is5 = 2A^,,. .... 

Quantum mechanics can violate inequality ([3]), but in 
order to do so, observables Oi and the state considered 
should be judiciously chosen. In order to check whether 
an inequality of the type of ^ allows for a non-locality 
test, one should maximize the left side of Q for an arbi- 
trary bipartite quantum state. If the maximum obtained 
value is greater than S, all states violating ([3]) are non- 
local. 

Before proceeding on testing the power of inequali- 
ties of the type of ^ for non-locality tests using the 
proposed molecular orientation-related observables, we 
demonstrate some simple general properties of ^ which 
are independent of the observables Oi. Such properties 
are useful since they help to simplify the numerical op- 
timization while giving some physical insight. We show 
that, thanks to these properties, the number of degrees of 



freedom to be considered in order to evaluate the max- 
imal value of the violation is decreased. We start by 
decomposing operators Oi (</>,) in terms of unitary trans- 
formations. Using the group property C/(0, + 0j) = 
U{(j)i)U{(f)j) we can write: 

(C/i(0i)(^J72((A2))(Oi(O)(8)O2(O))+ (4) 

Ol(0)«)O2((/)'2-(/'2) + 
Ol((/.'i-(/.l)®O2(0)- 

Oli^[ - C^l) ® 02(0^ - 02)(C/|(0l) ® C/l(02)), 

or, more succinctly, 

B^ {Ui{4>i)®U2{(j>2)) (5) 

S(0, 0, 01 - , </>2 - cl>'2)iUl (01 ) ® Ul (02)) , 

with 

S(0, 0, 01 - 0'i, 02 - 0^) = Oi(0) ® 02(0) + (6) 

Oi(o)(» 02(0^-02) + 

Oi(0'i-0i)®O2(O)- 

Ol(0'i~0l)®O2(0^-02), 

The interest of such decomposition is that operator 
S(0,0, 01 — 01,02 — 02) depends only on two variables 
and possesses the same spectrum as B. Thus, S(0, 0, 0i — 
01, 02 — 02) a-nd B share the same entanglement witness- 
ing properties. Of course, B and ,6(0, 0, 0i — 0'i, 02 — 02) 
do not detect the same entangled states. Nevertheless, 
they are connected by a local unitary transformation. 

The first question we wish to answer concerns the use- 
fulness of orientation correlation measurement for non- 
locality tests. In order to answer that, it is enough to 
determine the norm of B{Q, 0, 0i — 0i, 02 — 02)- 

In the next section, we will study a specific example of 
an operator Oi. 2(01,2), and the physical meaning of the 
previous results will be made explicit. 

III. ORIENTATION BASED NON-LOCALITY 
TEST 

A. Molecular orientation and correlation 

We introduce now an essential ingredient for the in- 
equalities studied in the present paper, which is molec- 
ular spatial orientation. Molecular orientation is an ob- 
servable, and its classical correspondent is the spatial ori- 
entation of the molecular inter-atomic axis with respect 
to some reference frame. By supposing that the molecules 
are addressed and manipulated by a linearly polarized 
laser field, we can define the laser's polarization axis as 
z, and take it as a reference for molecular orientation. 
In this case, the orientation of a given molecular state 
can be defined as the average value of operator cos^, 
where is the angle relative to the z axis. Note that 
here, cos(0) is taken as an operator, in perfect analogy 
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to the position operator and related functions. Orienta- 
tion as defined above can be experimentally measured, 
as discussed below. 

Since we are interested on probing properties related 
to a two party system, our system is composed by two 
molecules. We suppose that they behave like rigid ro- 
tors that can freely evolve. Their state depends thus on 
the time t. The Hamiltonian describing each individual 
molecule's free evolution is Hi ~ Jf/fi^, where i ~ 1,2. 
It is expressed in units of the rotational energy. J,; is 
the angular momentum operator and therefore the as- 
sociated evolution operator is given by Ui{t) ~ g-^Hit^ 
where time t is written in units of the rotational period. 
Ui{t) is therefore time-periodic with period 1. 

For each molecule of the bipartite set, the orientation 
at time t is defined as the expectation value of the Oi (t) = 
U2^{t)'S>U^^{t)cos{0^)Ui{t)'»U2it) operator, {di{t)) = 

(■0010^(^)1-00), where |?/'o) is the initial state of the system. 
Orientation correlations between two particles are given 
by the average value of {C{ti,t2)} = {Oi{ti) 02(^2)), 
and this quantity can be measured at different times ti 
and t2 for each molecule. Operator cos^ is useful for 
entanglement detection since it "mixes" different values 
of j, without affecting their projection m. Previous works 
have considered correlations between different values of 
the projection m of a given (fixed) value of j [Til. [3l|. 

With an arbitrary accuracy, each molecule's state (sub- 
scripts have been omitted) can be considered to reside in 
a finite dimensional Hilbert space TC generated by the 
basis set m); < j < jmax, \m\ < j}, where |j, to) are 
the eigenstates of and J^. Note that TC has dimen- 
sion (jniax + 1)^. The corresponding wavefunctions are 
{9,ip\j,m) = Yjm{0, ^p), the spherical harmonics. In the 
finite space 7^^^'"='^-', the cos 9 operator is characterized 
by a discrete, non degenerate spectrum of eigenvalues 
A„, with corresponding eigenvectors |A„), also called ori- 
entation eigenstates. The two maximally oriented states 
|-|-) and |— ) are the two eigenstates corresponding to the 
extreme eigenvalues zLX^ax, where Xmax = Max„(A„). 
In the particular case of jmax = 1 and to = 0, max- 
imally oriented states can be written in the basis of 
the angular momentum eigenstates |0, 0) and |1,0) as 
1+) = V^(|l,0) + |0,0))and|-) = yT72(|l,0)-|0,0)). 
In this particular case, cos 9 = y/l/ Sax and we recover 
the results of a two-level system. In this context, the 
free evolution operator [/,; (i) changes the orientation of a 
state, since orientation eigenstates arc not eigenstates of 
the free Hamiltonian. Time evolution creates superpo- 
sitions of orientation eigenstates, exactly as it happens 
when one projects photon's polarization with polarizers. 



B. The inequalities 

We can now combine all the ingredients to build the 
orientation based Bell-type inequalities. We start with a 
system composed by two molecules, and measurements 




FIG. 1: (Color online) Maximal value of (Bi) as a function of 
ti and t2 in units of the rotational period. Left jz-axis: highest 
eigenvalue /3i. Right z-axis: relative violation bi defined by 
Eq. Uni). (a): jmax = 1, ?n = 0, (b): jmax = 5, ?n = 0. 



independently performed in each one of them should be 
combined in order to tell whether the total molecular 
state violates or not local realism. The two molecule 
bipartite state \ipo) has been created at a given time 
to after which it freely evolves. We suppose that, af- 
ter to, there is no interaction between molecules. The 
total molecular state is thus given by the wavefunction 
ip{9i,92,(pi,(p2,t) = {9i,92,(pi,(p2\ipit)), where 9^ and (p^ 
denote here the polar and azimuthal spherical coordi- 
nates in the laboratory frame. 

This state is the one whose non-local properties are 
to be checked. We can use it to compute the average 
values of {C{ti,t2)) and combine such correlations taken 
at different times in a way analog to Eq.( [T|) and 
defining the operator 

6i(ti,t2,t'i,4)= (7) 
C{t^,t2)+C{t,,t'^) + C{t[,t2)-C{t[,t'^). 

For a local theory (LT), it obeys an inequality similar to 
Eq. ©: 



|(6i)lt| <2(A™a,)2; V(t„tOe 



(8) 



Without loss of generality, we have assumed that each 
particle state resides in the same finite dimensional 
Hilbert space 

7^(imax)^ Notice that, in Equation jZl), time 
plays the role of polarisers in Bell inequalities based in 
the photonic polarisation. Other CHSH inequalities us- 
ing the time evolution instead of polarizers were studied 
in the literature in very different contexts: in [26l. [3^. Issj 
they allow the detection of entanglement between prod- 
ucts of decaying mesons. In [3^.[35j. they reveal quantum 
properties of single particles. 
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FIG. 2: Relative violation of the inequality (|8]) as a function 
of the dimension of the subspace for m = 0. Blue: relative 
violation with respect to the locality threshold when only the 
finite dimension subspace is considered. Red: relative viola- 
tion when the infinite dimensional threshold of 2 is considered. 



We note that ([5]) is valid for all possible values of 
jmax, and that it can, in particular, be extended to the 
limit j'niax +00, in which case the spectrum of cos 9 
forms a continuum. An interesting characteristic of the 
separability threshold ([5]) is its dependence on Xmax- 

We will see in the next subsection how the general 
properties of CHSH-type inequalities can be used in the 
specific case discussed in this paper to show that the 
studied inequalities can be violated. 



C. Reference Frame and Temporal Origin 

A straightforward application of the results derived in 
section Hill consists of showing that, in order to study the 
spectrum of Bi , we can, without loss of generality, numer- 
ically diagonalize it in the specific case of t[ = t'2 = 0. 
This result is due to the fact that both operators are 
related by a local unitary transformation. Using the no- 
tation of Section II, we have here that ti = 0^ — (pA 
and t2 = 4>'g — ipB- Another interesting application of 
the general results of Section II consists of showing that 
the inequalities discussed here allow for local realism vi- 
olation even in the case where observers A and B have 
different time origins. The same happens for the spatial 
reference frame: violation is independent of any previous 
agreement between observers. However, different tempo- 
ral origins and reference frames correspond to different 
non-local detected states, that are related to each other 
by local unitary transformations. 

We start by discussing in more details the time origin 
chosen by both observers. Usually, the orientation Bell- 
type inequalities depend on four times of measurements, 
as defined in Eq. ([7]). However, as pointed out previously, 
different times correspond to the application of different 
unitary transformations. We can thus apply the results 



of section IIIIl identifying the general operator O to the 
operator cos6'. This leads to the inequalities: 

{BiihMAA)) = (Si(0,0,f'i-ti,4-<2)) <2AL.,, 

(9) 

The first identity shows that states maximally violat- 
ing the Bell- type inequality for ti,t2,t'i,t'2, defined as 
\smaxiti,t2,t[,t'2)) Can be obtained by the one maxi- 
mally violating it for 0, 0, t'j^ — ii, — that we will call 
|smai:(0, 0, t']^ — ti, ^2 — ^2)) by the application of the trans- 
formation U{ti)'» U{t2)\smax{0, 0, t[ - <i, <2 - ^2))- Tem- 
poral uncertainties of ri , T2 for each one of the molecules 
can always be translated as the application of the uni- 
tary operator U{ti) (E) U(t2), so that their only effect is 
to change the eigenstatcs of operators Bi by the same 
transformation. Violation can thus still be observed, and 
the subspace of states violating local realism are obtained 
by a simple unitary transformation on the original sub- 
space. 

The same type of argument can be used for uncertain- 
ties of the reference frame for each observer. All reference 
frames arc connected by local unitary transformations 
describing rotations about some direction of space, and 
results for different references frames are connected by 
these same unitary transformations. 

We apply these results to simplify the investigation of 
possible violations of inequality ([8]) . 



D. Results 

For a given value of jmax, we have numerically diag- 
onalized operator Bi{0,0,ti,t2), and obtained for each 
(ti,i2), its highest eigenvalue Pi{ti,t2), which gives the 
maximal value of (Bi) (maximal violation of ([5])). This 
quantity depends on the dimensionality of the system, 
and we compare the amplitude of the violation when dif- 
ferent values of jmax are used by defining the relative 
violation 



6l(tl,t2) 



max J 

2 ( Amaa; )^ 



(10) 



Results for different values of jmax and m = are shown 
in Figure [TJ We can see that the proposed inequahties 
are violated for a significant region of parameters <i and 
t2- Figure [T] also calls one's attention because of its sym- 
metries. The central symmetry with respect to the point 
ti = t2 = 0.5 corresponds to the time reversal symme- 
try. One can also easily understand the mirror symmetry 
about the ti = t2 line with the help of the particle ex- 
change symmetry of operator Bi . 

We have shown in [22| that Eq. ([5]) can be violated 
by a number of pure states, and numerically studied the 
relative violation (jlOp with increasing dimension. The re- 
sults, shown in Figure ([2]) were obtained in a particular 
case, where both molecules had a vanishing angular mo- 
mentum z axis projection (m = for both molecules). 
The effects of considering different values of m will be 
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discussed in the following. We focus here in the to = 
case in order to discuss the behaviour of violation of ([5]) 
with respect to the dimensionality of the system. We see 
that the maximal relative violation 6i(ti,<2) decreases 
with the dimensionality of the system for low dimensions 
and then starts smoothly growing starting from jmax = 4 
(blue bars, Figure Up to now, wc have not found 

an asymtoptical numerical behaviour for jmax oo. We 
have shown that the maximal violation is bounded by 3. 
However, it is still unknown if this value can be reached. 
Since the exact behaviour of the violation of inequality 
([5]) in infinite dimension is still not known, one perti- 
nent question is whether the proposed inequality is still 
violated in the limit jmax oo for physically accept- 
able states. We can see that it is indeed the case by 
calculating the violation relative to the classical bound 
obtained for an infinite dimensional system. In this case, 
the classical threshold separating local theory from non- 
local ones is 2, since Xmax = 1- Figure shows the 
results of such violation (pink bars). We see that states 
violating the infinite dimensional classical threshold can 
be found when one considers subspaces with dimensions 
higher than N = jmax -I- 1 > 5 for m = 0. In this sub- 
space, we can see that the violation relative to the infinite 
dimensional threshold is still small. However, by increas- 
ing the size of the subspace and going to jmax = 10, we 
can see that this relative violation increases, and it is not 
negligible when compared to the two dimensional case, 
for instance. Notice that the relative violation with re- 
spect to the infinite dimensional subspace and the one rel- 
ative to the restricted subspace approach with increasing 
dimension, as one should expect. Violating the inequal- 
ities (|S]) in the infinite dimensional limit with entangled 
states of low dimension is a surprising result that proves 
that ^ can be violated for all possible values of Xmax- 

A natural question is what states maximally vio- 
late (|S|). Since we are dealing here with correlations 
in orientation, one could expect that maximally ori- 
ented entangled states of the form -^{\Xmax, Xmax) ± 
\—Xmaxi —Xmax)) SLTC tliosc that maximally violate them. 
However, this is not the case, and such state, except for 
the case jmax = 1 and to = do not violate ([5]) at all. 
Nevertheless, states maximally violating ([8]) have most of 
their population in highly oriented states. For jmax = 5, 
states \Xmax,Xmax) and \-Xmax,-Xmax) Carry, each one 
of them, 36% of the population. In addition, states max- 
imally violating ([8]) are not maximally entangled (except 
in the case jmax = 1, a result which is not surprising and 
has been observed for other Bell-type inequalities, as in 
??, for example. Using the maximally violating state, 
we have computed the reduced density matrix entropy 
S = —Tr[pilog Pi], where pi is the reduced density ma- 
trix with respect to one of the two entangled molecules. 
In order to compare S for different values of jmax, we 
have normalised it with respect to log {jmax + 1), which 
is the maximal value of the entropy in a subspace of di- 
mension Jmax + 1- The reduced density matrix entropy 
is a measure of entanglement for pure states, as the ones 
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FIG. 3: Normalised entropy of the maximally violating state 
as a function of the subspace dimensionalty jmax (dots). Con- 
tinuous (green) line represents log 2/ log (jmax + !)• 

we are considering here. The results are shown in Fig- 
ured!]), together with a plot of log 2/ log (jmax + !)■ We 
see that the two plots almost completely match, showing 
that maximally violating states are very close to entan- 
gled states involving only two orthogonal states of each 
molecule. We have checked that the two states involved 
are the maximally oriented states. 

The results presented above were obtained for the spe- 
cific case of TO = for both molecules. Allowing to to 
take a value different from to = does not bring quali- 
tative changes to our results: the classical threshold will 
still be given by the case m ~ 0, since it corresponds to 
the maximal value of (Bi) for a local state. However, 
the dimension of the subspace considered depends on to. 
For a given fixed to, it is given by jmax — |™| + 1- The 
value of the maximal orientation is thus determined not 
only by the dimension of the subspace, but also on the 
specific value of m. As a consequence, the maximal value 
of (Bi) is determined by the same parameters. For in- 
stance, if the considered state is a superposition state 
of different to's, each one of this subspaces will lead to 
different contributions when computing (Bi). In partic- 
ular, those for which jmax ~ |™| give a null contribution, 
and the maximal possible contribution decreases as |to| 
increases for each molecule. Physically, this is related to 
the fact that states with high values of the angular mo- 
mentum projection |to| arc less oriented than those for 
which J ;g> |m|. As a conclusion, by considering different 
to's, it is still possible to violate ([8]), even if fulfilling the 
required conditions for it becomes harder. 

E. Non-pure states and the effect of noise 

We now formalize the conditions a non-pure state 
should satisfy in order to violate Eq. ©. This is useful 
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to estimate the effect of noise in our system. A non-pure 
bipartite state is described by the density matrix pi_2i 
where the subscripts 1 and 2 refer to each one of the 
observers. The average values of ([7]) can be obtained 
from pi^2 by using (C(ii,f2)) = (Oi(ti) ® 02(i2)> = 
Tr[pi^20i(ii) ® 02(^2)]- Defining the Bi eigenstates, by 
{|si)}, with i ranging form to (jmaa: + l)^ if one assumes 
that each particle's subspace have the same dimension 
and m = 0, it is clear that the density matrix pi^2 can be 
expressed in this basis. By doing so, we have that 

Tr[Sipi,2] =^P^s„ (11) 

i 

where si is the eigenvalue associated to the eigenstate \si) 
of Bi.pi = {si\pi^2\si} the statistical weight of each one of 
such eigenstates and N = jmax + 1 is the dimensionality 
of each molecule's subspace. This means, that for a local 
realistic theory, one should obey 

Y,P^s, < 2AL.- (12) 

i=l 

An example of non-pure state is given by 

Pi, 2 = + (1 - PN)\Sm.ax){Sniax\, (13) 

where 1 is the iV^ x identity matrix, \smax) is the 
eigenstate of Bi with maximal eigenvalue Smax and Pn 
is the probability of the state to be in a complete mixture. 
This type of state has been studied in (2^ and can illus- 
trate the presence of noise in the preparation of a state 
|s?nai:)niaximally violating a Bell- type inequality. They 
defined that the robustness of a Bell-type test with re- 
spect to noise is measured by the maximal allowed value 
of Pn still leading to locality violation. In the notation 

of Eq. (ini), {s^ax\p\Smax) = Pmax = Pn/N^ + (1 " Pn) 

and (sj ^ s„iax\p\si 7^ Smax) = Pi = Pn/N"^- By cal- 
culating the expectation value of Bi using one gets 
that local realism should obey: 

E ]^ + (1 - PN)Smax < 2X1,,. (14) 
1=1 

Since Bi is traceless, — 0, so the inequality is 

violated for Pjv < 1 — ^X^^^/smax- A plot of Pn as a 
function of the considered subspace is given in Figure ([4]) . 
Notice that, when a two-level system is considered, the 
known result of P2 = 1 ^ is recovered. Contrary to 
other Bell-type inequalities for which Pn was optimised 
[2^, the maximal allowed value of noise in our system 
decreases with dimensionality up to jmax = 5 and then 
starts increasing again, but doesn't change significantly 
in the range of dimensions that we have calculated. This 
result shows that the non-local properties of the max- 
imally violating state are quite robust with respect to 
noise. 
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F. Dichotomisation 

We have shown that the operator defined by Eq. ^ 
allows not only for the realization of BcU-type tests in 
finite angular momentum subspaccs. but also when the 
size of the subspace is not a priori known. This proves 
that, even in the case where we consider an infinite di- 
mensional space, our inequalities can be violated for some 
states. However, as shown previously, the contrast of the 
maximal violation depends on the value of jmax ■ With 
dichotomizing procedure, a high dimensional system is 
transformed into an effective two level one by splitting in 
two the set of measurement outcomes, and the contrast 
is kept constant, equal to the one for a two- level system, 
irrespectively of jmax- In our case, we can dichotomise 
as follows: states \(p), for which (cos6')<p > 0, arc said to 
be positively oriented, while those for which {cos9)ip < 
are considered as negatively oriented. Wc define the as- 
sociated projectors: U± ~ J2x±\^±){^±\ which project 
states in the subspace of positive (negative) orientation. 
A± are the positive (negative) cos 9 eigenvalues in the 
given finite dimensional subspace and |A±) are the corre- 
sponding eigenstates. The measured observable for each 
molecule i is then Hi = 11+ — n_ . For a given single- 
molecule state \ip) = Y,x+ CA+|A+) +J2x_ ca_|A_), (11)^ 
can take any value in the interval [-f 1, — 1]. The total two 
molecule observable is 11 = Hi (X) 112. We refer, as before, 
to two molecule correlation measurements realized at two 
different times, using n(ii,i2) = ^i{ti) (Sill2{t2) where 
Ili(ti) — U^^{ti)UiUi{ti). In analogy with Eq.[71 we now 
define the operator 

S2 = n(ti,i2) + n(ti,i;,) + n(t;,t2)-n(f;,t'2). (i5) 

Since n(ti, tj)'^ = 1, one can show that the highest value 
(B2) can reach is given by the Cirel'son bound 2^/2 [29l |. 
Also, for a local theory, we have 

|(S2)lt|<2; V(ti,t2)eM2. (16) 
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FIG. 5: (Color online) Same as Fig. [Ub, but for {B2) 



As in the case of Bi, we can define the relative viola- 
tion b2(ti,t2) as a function of the maximal eigenvalue 
/32(<i,i2) of i32(ii,i2) (see Equation (dHD). The interest 
of dichotomisation is that the maximal value of the rel- 
ative violation is always given by &2(ii,i2) = — 1, 
independently of the dimension of the subspace. Fig. O 
shows the maximum value /32(iiji2) of {B2{ti,t2)) as a 
function of ti and t2 for jmax = 5, while for jmax = 1 we 
obtain trivially the same result as with 6i(ti, (Fig.[TJa 
right z-axis). We see in Fig. Othat the dichotomisation 
procedure not only keeps the contrast of the maximal vi- 
olation constant, but also allows for high violation for a 
wider range of values of and ^2- 



IV. DISCUSSION ON POSSIBLE 
EXPERIMENTAL IMPLEMENTATIONS 

Orientation entanglement between two molecules can 
be created, for instance, with the help of the dipolar in- 
teraction [ij. A possible experimental scenario consists 
of two trapped molecules that can be submitted to spa- 
tial displacements, as it is currently done with atoms 
[2^. By putting molecules close enough to each other, 
and with the help of laser manipulation, one can tai- 
lor rotational entangled states. In order to realise Bell- 
type tests, molecules would have to be separated and 
measurements performed independently in each one of 
them. The proposed scenario is realistic and there has 
been rapid progress on cooling, trapping and manipu- 
lating polar molecules, especially those formed with two 
alcali atoms. We believe that such manipulations are go- 
ing to be shown in laboratories in a near future. With 
respect to the tests wc wish to perform in the present 
paper, there are some crucial points to analyse for the 
experimental implementation of the proposed Bell-type 
inequality test. One of them concerns trapping condi- 
tions. Traps should allow molecules to make spatial dis- 
placements so that they are sufficiently far apart and so 
that signaling of one observer's results to the other is 
not possible during at least one molecular rotational pe- 
riod. This condition is attained for an inter-molecular 
distance L > cT, where c is the speed of light and T 
is the molecular rotational period. To put some num- 



bers on it, the molecular period is usually of the order of 
10~^^s, giving a lower bound to L k, 10~^m. Even if at 
this distance dipolar interaction between molecules is not 
completely negligible, measurements on both molecules 
are performed during the same rotational period, which 
is much faster than the characteristic interaction time 
related to dipolar coupling at this distance (see for 
example, for a detailed discussion on the relevant time 
scales) . Optical traps seem promising candidates to trap 
and to displace molecules. Optical tweezers [11], opti- 
cal lattices and optical conveyor belts [2^ are some ex- 
amples of dipole force based optical traps allowing for 
atomic trapping and coherence preserving displacements 
m, m, [13] • In the case of optical tweezers, highly fo- 
cused beams are used to trap individual atoms [36| . and 
we can, in principle^Joad two atoms, each one of them in 
an individual trap [38|. In optical conveyor belts, a sta- 
tionary field made of two counter propagating beams is 
used to trap atoms, that can be displaced over distances 
of « lOnm by changing the relative detuning between 
the laser beams [2^. Going from atoms to molecules in 
this type of experimental system could in principle be 
done by the usual photo association techniques, already 
demonstrated in the context of optical lattices [l^ , where 
diatomic molecules of two alcali atoms of the same specie 
have been produced. Once molecules are created and 
trapped, they interact with each other by dipole interac- 
tion and a non-local state may be created, as discussed 
in the previous section. Molecules are thus separated un- 
til dipole interaction is negligible, and measurements can 
be performed. One natural question concerns the effect 
of the trapping potential itself to the orientation of the 
molecular state. Optical traps are based on an applied 
non resonant electric field that interacts with rotational 
levels. We suppose that molecules in our setup are cre- 
ated in their ground electronic level. The optical trap is 
created by coupling non-resonantly this electronic level 
to the first excited one. The detuning between laser and 
the electronic transition is 5. In usual experiments on op- 
tical tweezers, for instance,, 5 can vary vary in the range 
~ 10 — lO^GHz, providing photon emission rates from the 
excited electronic state in the range of 0.1 — 100 MHz. 
These figures are of the same order of magnitude of co- 
herence lifetime for rotational levels, while the range of 
variation of 5 is also of the same order of magnitude of 
the frequency difference between neighboring rotational 
levels. In order to estimate various effects of trapping 
lasers, we will put ourselves in a close to realistic con- 
figuration in which the rotational frequency is « lOGHz 
and 5 « lO^GHz. Each rotational level has a different 
energy, making the effective detuning j dependent. This 
corresponds to adding a factor 5j ~ Bj{j + l)/h to the 
electronic detuning S, taken with respect to the rota- 
tional ground state. Because each rotational level has a 
detuning which depends on the value of j, there will be a 
phase difference between each rotational level due to the 
non resonant coupling to light. In the case of a harmonic 
trapping potential, this phase is given by e~*^^^ 



9 



if only one electronic transition is considered. This cor- 
responds to a local action on each molecule that does 
not play a significant role on non-locality tests, since it 
maps one state to another one with the same degree of 
entanglement. As shown before, by properly choosing 
the time where measurements are performed, such effect 
can be compensated. This means that the parameters for 
which maximal violation occurs may be modified. As it 
can be shown, if a given state violates the inequalities of 
the type of ^ and pB)) . a collection of states connected 
to it by local unitary transformations will also violate the 
same inequalities, even if for different parameters. Also, 
it is important to notice that such dephasing effects can 
be completely neglected depending on the precise circum- 
stances under which the experiment is performed. This 
happens because dephasing occurs on a time scale much 
longer than the rotational period, so that if measure- 
ments are performed in this interval, it will not affect the 
expected results. 

Finally, we address the question of orientation mea- 
surement. A possibility is to use Coulomb explo- 
sion, a destructive techniques alread y e mployed for 
molecular orientation measurements [30|. After a 
quasi-instantaneous (when compared to the molecu- 
lar rotational period), laser induced dissociative multi- 
ionisation, the molecular fragments are recorded at dif- 
ferent directions of space, in the case of the orientation 
Bell- type inequalities, or at different hemispheres, in the 
case of the dichotomized inequalities. This technique is 
still experimentally challenging since it demands highly 
efficient single atom detectors. Usually this technique is 
employed for molecular ensembles, and the detection ef- 
ficiency is less determinant than in our case. Another 
possibility is to detect the orientation of single molecules 
optically, as realized in [2l| . Fluorescence intensity when 
a single molecule is excited by a linearly polarized laser 
beam that can change polarization in time, can reveal 
molecular orientation since the scalar product between 
the field polarization and the molecular dipolc depends 
on its orientation. The phase of the fluorescence inten- 
sity with respect to the exciting laser beam when its po- 
larization is turned in time, depends on the molecular 
orientation, as observed in [2]| . 



V. CONCLUSION 

We have extensively discussed some important prop- 
erties of a recently proposed Bell-type inequality based 
on the measurement of a continuous, bounded, observ- 
able, which is molecular orientation. We have discussed 
some important properties of it, as its symmetries and 
its non-dependency on a specific choice of common tem- 
poral origin or reference frame. The role of noise in our 
system, that would transform pure states into statistical 
mixtures was also analysed, as well as how violation de- 
pends on the amount of noise. We have also discussed 
some important conditions that an eventual experimen- 
tal set-up should satisfy and the influence of physical pa- 
rameters as light forces and residual interaction between 
molecules on the proposed measurements. It seems, from 
our analysis, that once the necessary degree of advance- 
ment is attained by experimental set-ups, our proposal 
is realistic and should not present major difficulties. Our 
results open the perspective of entanglement detection 
and non-locality tests for high angular momentum sys- 
tems in atomic and molecular physics. 

One interesting aspect of the derived inequalities is the 
fact that they assume the simple form of CHSH inequali- 
ties and deal with continuous variables at the same time. 
This happens because the observable that is measured for 
each particle has a bounded spectrum, naturally limiting 
the locality threshold and the norm of the Bell operator. 
Notice that, in spite of having discussed here the spe- 
cific case of molecular orientation, the same type of in- 
equalities could be derived for other continuous bounded 
observables, easily measurable in other physical contexts. 
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